Block Diagonalization using SRG Flow Equations 
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By choosing appropriate generators for the Similarity Renormahzation Group (SRG) flow equa- 
tions, different patterns of decoupling in a Hamiltonian can be achieved. Sharp and smooth block- 
diagonal forms of phase-shift equivalent nucleon-nucleon potentials in momentum space are gener- 
ated as examples and compared to analogous low- momentum interactions ("Mowfc")- 
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The Similarity Renormalization Group (SRG) [U |2l 
[3] applied to inter-nucleon interactions is a continuous 
series of unitary transformations implemented as a flow 
equation for the evolving Hamiltonian Hs^ 

^ = [r?„F,] = [[G„F,],F,]. (1) 

Here 5 is a flow parameter and the flow operator Gg spec- 
ifies the type of SRG [4 . Decoupling between low-energy 
and high-energy matrix elements is naturally achieved in 
a momentum basis by choosing a momentum-diagonal 
flow operator such as the kinetic energy T^ei or the diag- 
onal of Hs] either drives the Hamiltonian toward hand- 
diagonal form. This decoupling leads to dramatically im- 
proved variational convergence in few-body nuclear sys- 
tems compared to unevolved phenomenological or chiral 
EFT potentials [110. 

Renormalization Group (RG) methods that evolve NN 
interactions with a sharp or smooth cutoff in relative 
momentum, known generically as Mowfe, rely on the in- 
variance of the two-nucleon T matrix [3 [8] . These ap- 
proaches achieve a block- diagonal form characterized by 
a cutoff A (see left plots in Figs, [l] and [2|. As usually im- 
plemented they set the high-momentum matrix elements 
to zero but this is not required. 

Block-diagonal decoupling of the sharp Fiow k form can 
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FIG. 1: (Color online) Comparison of momentum-space Mowfe 
(left) and SRG (right) block- diagonal potentials with A = 
2fm~^ evolved from an N^LO ^Si potential 11 . The color 
axis is in fm. 



be generated using SRG flow equations by choosing a 
block-diagonal flow operator [9l [10] , 

r — ( P^sP \ _ „bd (r.\ 

QHsQ )=^s . (2) 

with projection operators P and Q = 1 — P. In a partial- 
wave momentum representation, P and Q are step func- 
tions defined by a sharp cutoff A on relative momenta. 
This choice for Gg, which means that 77^ is non-zero only 
where Gg is zero, suppresses off-diagonal matrix elements 
such that the Hamiltonian approaches a block-diagonal 
form as s increases. If one considers a measure of the 
off-diagonal coupling of the Hamiltonian, 

TT[{QHsPy{QHsP)] = Tt[PHsQHsP] ^ , (3) 

then its derivative is easily evaluated by applying the 
SRG equation, Eq. ([T]): 

^Tt[PHsQHsP] 
as 

= TT[Pr]sQ{QHsQHsP - QHsPH.P)] 

+ Tt[{PHsPHsQ - PHsQHsQ)Qr]sP] 
= -2Tr[(Q77,P)t(Q7^,P)] ^0. (4) 

Thus, the off-diagonal QHqP block will decrease in gen- 
eral as s increases [9l [10] . 




FIG. 2: (Color online) Comparison of momentum-space Mowfe 
(left) and SRG (right) block-diagonal potentials with A = 
2fm~^ evolved from an N^LO ^Si potential [iT]- The color 
and z axes are in fm. 
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The right plots in Figs. [T] and [2] result from evolving the 
N^LO potential from Ref. [11 using the block-diagonal 
Gs of Eq. ^ with A = 2fm"^ until A = l/s^^^ = 
0.5 fm~^. The agreement between Viowk and SRG poten- 
tials for momenta below A is striking. A similar degree of 
universality is found in the other partial waves. Deriving 
an explicit connection between these approaches is the 
topic of an ongoing investigation. 

The evolution with A of two representative partial 
waves (^Si and ^Pi) are shown in Figs, [s] and H The 
evolution of the "off-diagonal" matrix elements (mean- 
ing those outside the PHgP and QHgQ blocks) can be 
roughly understood from the dominance of the kinetic 
energy on the diagonal. Let the indices p and q run 
over indices of the momentum states in the P and Q 
spaces, respectively. To good approximation we can re- 
place PHgP and QHsQ by their eigenvalues Ep and Eg 
in the SRG equations, yielding [9l fTO] 

-^hpq ~ TjpqEq — EpTjpq = —{Ep — Eq) Tjpq (5) 

and 

Tjpq Ephpq hpqEq = {Ep Eq) Jlpq . (6) 

Combining these two results, we have the evolution of 
any off-diagonal matrix element: 

-^hpq ^ -{Ep - Eq)'^ hpq . (7) 

In the NN case we can replace the eigenvalues by those 



for the relative kinetic energy, giving an explicit solution 
/ij,,(s)«/ip,(0)e-«(^''-^^)' (8) 

with Cp = /M. Thus the off-diagonal elements go to 
zero with the energy differences just like with the SRG 
with Trei; one can see the width of order l/y/s = A^ in 
the /c^ plots of the evolving potential in Figs, [s] and |4j 

While in principle the evolution to a sharp block- 
diagonal form means going to 5 = oo (A = 0), in practice 
we need only take s as large as needed to quantitatively 
achieve the decoupling implied by Eq. ([8|. Furthermore, 
it should hold for more general definitions of P and Q. 
To smooth out the cutoff, we can introduce a smooth 
regulator /a, which we take here to be an exponential 
form: 

fA{k) = e-f'^V^^)" , (9) 

with n an integer. For Mow/c potentials, typical values 
used are n = 4 and n = 8 (the latter is considerably 
sharper but still numerically robust). By replacing H^^ 
with 

Gs = fAHjA + {l-h)H,{l-fA), (10) 

we get a smooth block-diagonal potential. 

A representative example with A = 2 fm~^ and n = 4 is 
shown in Fig.js] We can evolve to A = 1.5 fm~^ without a 
problem. For smaller A the overlap of the P and Q spaces 
becomes significant and the potential becomes distorted. 
This distortion indicates that there is no further benefit 
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FIG. 5: (Color online) Evolution of the ^Si partial wave with a smooth (n = 4) block- diagonal flow equation with A = 2.0 fm~^, 
starting with the N^LO potential from Ref. [11^. The flow parameter A is 3, 2, 1.5, and lfm~^. The axes are in units of 
from 0-11 fm~^. The color scale ranges from —0.5 to +0.5 fm as in Fig. [l] 



to evolving in A very far below A; in fact the decoupling 
worsens for A < A with a smooth regulator. 

Another type of SRG that is second-order exact and 
yields similar block diagonalization is defined by 



Vs 



[T.PVsQ^QVsP] 



(11) 



which can be implemented with P ^ /a and Q 
(1— /a), with /a either sharp or smooth. We can also con- 
sider bizarre choices for /a in Eq. (10), such as defining 
it to be zero out to Aiower, then unity out to A, and then 
zero above that. This means that 1 — /a defines both low 
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FIG. 6: (Color online) Evolved SRG potentials starting from 
Argonne vig in the ^So and ^Pi partial waves to A = lfm~^ 
using a bizarre choice for Gs (see text). The color and z axes 
are in fm. 



and high-momentum blocks and the region that is driven 
to zero consists of several rectangles. Results for two par- 
tial waves starting from the Argonne vis potential [12] 
are shown in Fig. [6| Despite the strange appearence, 
these remain unitary transformations of the original po- 
tential, with phase shifts and other NN observables the 
same as with the original potential. These choices pro- 
vide a proof-of-principle that the decoupled regions can 
be tailored to the physics problem at hand. 
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FIG. 7: (Color online) Phase shifts for the ^Si partial wave 
from an initial N'^LO potential and the evolved sharp SRG 
block- diagonal potential with A = 2fm~^ at various A, in 
each case with the potential set identically to zero above A. 

Definitive tests of decoupling for NN observables are 
now possible for Flow k potentials since the unitary trans- 
formation of the SRG guarantees that no physics is 
lost. For example, in Figs, [t] and [s] we show ^Si phase 
shifts from an SRG sharp block diagonalization with 
A = 2fm~^ for two different potentials. The phase shifts 
are calculated with the potentials cut sharply at A. That 
is, the matrix elements of the potential are set to zero 
above that point. The improved decoupling as A de- 
creases is evident in each case. By A = lfm~^ in Fig. [7| 
the unevolved and evolved curves are indistinguishable 
to the width of the line up to about 300 MeV. 
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FIG. 8: (Color online) Same as Fig.|7|but with Argonne vig 
as the initial potential [12]. 
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FIG. 9: (Color online) Errors in the phase shift at ^lab = 
100 MeV for the evolved sharp SRG block-diagonal potential 
with A = 2 fm~^ for a range of A's and a regulator with n = 8. 



In Fig.|9]we show a quantitative analysis of the decou- 
pling as in Ref. [13 . The figure shows the relative error 
of the phase shift at 100 MeV calculated with a poten- 
tial that is cut off by a smooth regulator as in Eq. (|9| 
at a series of values Acut • We observe the same universal 
decoupling behavior seen in Ref. [13 : a shoulder indicat- 
ing the perturbative decoupling region, where the slope 
matches the power 2n fixed by the smooth regulator. The 
onset of the shoulder in Acut decreases with A until it sat- 
urates for A somewhat below A, leaving the shoulder at 
Acut ^ A. Thus, as A ^ the decoupling scale is set by 
the cutoff A. 

In the more conventional SRG, where we use r]s = 
[T^Hs] = [T, Fg], it is easy to see that the evolution of 
the two-body potential in the two-particle system can 
be carried over directly to the three-particle system. In 
particular, it follows that the three-body potential does 
not depend on disconnected two-body parts [H If 
we could implement rjs as proposed here with analogous 
properties, we would have a tractable method for gener- 
ating Viowk three-body forces. While it seems possible to 
define Fock-space operators with projectors P and Q that 
will not have problems with disconnected parts, it is not 
yet clear whether full decoupling in the few-body space 
can be realized. Work on this problem is in progress. 
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